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Abstract. In this paper, subgraphs and complementary graphs are used to analyze the 
network synchronizability. Some sharp and attainable bounds are provided for the eigenratio of 
the network structural matrix, which characterizes the network synchronizability, especially when 
the network's corresponding graph has cycles, chains, bipartite graphs or product graphs as its 
subgraphs. 
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1 Introduction 

Synchronization of complex networks has attracted increasing attention from many scientists, for 
its important applications in various areas of physical and biological sciences [3"1 HI IS] \9\ [TOl [231 [26]. 
Synchronization is a ubiquitous phenomenon especially in social and biological networks where, 
quite often, it is desirable, e.g., in the consensus of multi-agent activities, while in some other cases 
it is undesirable, e.g., in traffic congestion [6], [TBI HH1 UHl E2], [231 E5] • The connection structure of a 
network plays an important role in its synchronization. Of particular interest in this concern is how 
the synchronizability depends on various structural parameters of the network, such as average dis- 
tance, clustering coefficient, coupling strength, degree distribution and weight distribution, among 
others. Some interesting results have been established for such important problems based on the 
notions of master stability function and synchronized region EJ [TTj [T2l fT3"j [TBI [2~T1 [29] . Some rela- 
tionships between synchronizability and structural parameters of complex networks have also been 
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reported, e.g., smaller average network distance does not necessarily mean better synchronizability 
|19j . therefore the betweenness centrality was proposed as a good indicator for synchronizability 
|10j . and two networks with the same degree sequence can have different synchronizabilities [27] . 
Some optimizing rules were established for the network synchronizability [6j. Moreover, how the 
network statistical properties influence the spectrum of the structural matrix was analyzed and re- 
ported in [2j. And the effects of graph operations on the network synchronizability were studied in 
PQ. Last but not least, complementary graphs were used to analyze the network synchronizability in 
[7]. Motivated by all these research works, this paper attempts to further establish some sharp up- 
per bounds for the network synchronizability based on the theory of subgraphs and complementary 
graphs. 

Consider a dynamical network consisting of N coupled identical nodes, with each node being an 
n-dimensional dynamical system, described by 

N 

Xi = f(xi) - c J ^2a ij H(x j ), i = 1,2, ••• ,N, (1) 

where X{ = (xn, Xi2, ■ ■ ■ ,Xi n ) E M n is the state vector of node i, /(•) : M n — » 1" is a smooth 
vector-valued function, constant c > represents the coupling strength, H(-) : M. n — > R n is called 
the inner linking function, and A = (aij)NxN is called the outer coupling matrix or structural 
matrix, which represents the coupling configuration of the entire network. This paper only considers 
the case that the network is diffusively connected, i.e., A is irreducible and its entries satisfy 
a u = ~~ Ylj=i j^i a iji * = 1,2,--- , N. Further, suppose that, if there is an edge between node i and 
node j, then = aj% = — 1, i.e., A is a Laplacian matrix corresponding to network (pQ). Therefore, 
is an eigenvalue of A with multiplicity 1, and all the other eigenvalues of A are strictly positive, 
which are denoted by 

= Al < A2 < A3 < • • • < Xn = Xmax- (2) 

The dynamical network (TjQ) is said to achieve (asymptotical) synchronization if x\(t) — > X2(t) — > 
■ ■ ■ — > xjv(t) - ► s(t), as t — > 00, where, because of the diffusive coupling configuration, the syn- 
chronous state s(t) € W 1 is a solution of an individual node, i.e., s(t) = f(s(t)). 

It is well known [3] that the eigenratio r(A) = j^- of network structural matrix A characterizes 
the network synchronizability: the larger the r(A), the better the synchronizability. In this paper, 
r(A) will be used as a synchronizability index for the networks with structural matrix A. The 
relationships between r(A) and network structural characteristics such as average distance, node 
betweenness, degree distribution, clustering coefficient, etc. have been studied [3 EJ [27l [28] . A 
well-known upper bound for the eigenratio r(A) in graph theory is r(A) < ^ min , where d m i n and 
dmax denote the smallest and largest degrees of the corresponding graph to the network [27J. To 
the best of our knowledge, besides this coarse bound, there are no sharp upper bounds available 



in the literature today. Some sharp upper bounds for r(A) are derived in this paper by using 
graph-theoretical methods. 

The rest of this paper is organized as follows. In Section 2, some necessary preliminaries in 
graph theory are provided. In Section 3, for a given graph G, some upper bounds are derived for 
the eigenratio r(G) when G has cycles as its subgraphs. In Section 4, some lower bounds for the 
largest eigenvalue of a graph are presented when G has other subgraphs, such as chains, bipartite 
graphs and product graphs. The paper is concluded by the last section. 

Throughout this paper, for any given undirected graph G, eigenvalues of G mean eigenvalues 
of its corresponding Laplacian matrix. Notations for graphs and their corresponding Laplacian 
matrices are not differentiated, and networks and their corresponding graphs are not distinguished, 
unless otherwise indicated. 



2 Preliminaries 

For a given graph G, let V and £ denote the sets of nodes and edges of G, respectively. A graph 
G\ is called an induced subgraph of G, if the node set Vi of G\ is a subset of V and the edges 
of G\ are all edges among nodes Vi in £ . The complementary graph of G, denoted by G c , is the 
graph containing all the nodes of G and all the edges that are not in G. In this paper, subgraphs 
and complementary graphs are used to discuss the network synchronizability. For this purpose, the 
following lemmas are needed. 

Lemma 1 [14, [15] For any given graph G of size N, its nonzero eigenvalues indexed as in ([2]) grow 
monotonically with the number of added edges; that is, for any added edge e, Xi(G + e) > Aj(G), 
i = l,... ,N. 

Lemma 2 |14j For any given connected graph G of size N, its largest eigenvalue An satisfies 
An > dmax + 1, with equality if and only if d max = N — 1, where d max is the maximum degree of 
G. 

Lemma 3 [8j [27] For any cycle Cn at N (> 4) nodes, its eigenvalues are given by /xi, • • • ,hn 
(not necessarily ordered) with [x\ = and 

sin(^) 

Mfe+i = 3 — . , Ls , k = 1, • • • , N - 1. 
sm(f ) 

By Lemma 3, one knows that, for any cycle Cn at N (> 4) nodes, if N is even, its largest 



sm ( 2N I 
Sm ( 2/V ) 



eigenvalue is 4; if N is odd, its largest eigenvalue is An(Cn) = 3 

Obviously, for cycles with odd numbers of nodes, the largest eigenvalue converges to 4 as the 
number of its lengths tends to +oo. For eigenvalues of graphs and complementary graphs, the 
following lemma is useful (see [14^ [To] and references therein) . 



Lemma 4 For any given graph G of size N, the following statements hold: 

(i) Aat(G), the largest eigenvalue of G, satisfies Xn(G) < N. 

(ii) \n(G) = N if and only if G c is disconnected. 

(iii) If G c is disconnected and has (exactly) q connected components, then the multiplicity of 
X N (G) =N is q- 1. 

(iv) Xi(G c ) + X N ^ +2 (G) =N, 2<i<N. 

For a given graph, generally its largest eigenvalue is easier to compute than the smallest nonzero 
eigenvalue. Hence, by Lemma 4, one can obtain the smallest nonzero eigenvalue by computing the 
largest eigenvalue of its complementary graph. 

Corollary 1 For any given graph G of size N, if its second smallest eigenvalue equals its smallest 
node degree, i.e., \2(G) = d m i n (G), then G or G c is disconnected; if A 2 (G) > d m i n (G), then G is a 
complete graph; if both G and G c are connected, then A 2 (G) < d m i n (G). 

Proof If G c is connected, then X max (G c ) = N — A 2 (G) = N - d mm (Gi) = d max (G c ) + 1. By 
Lemma 2, d max {G c ) = N — 1, so G is disconnected. Further, if A 2 (G) > d m i n (G), by Lemma 
4, Atv(G c ) < d maa; + 1. Combining with Lemma 2, G c can only have isolated nodes, i.e., G is a 
complete graph. The third statement holds obviously. □ 

Combining Lemmas 1 and 2 with Corollary 1, one can easily get the following result. 

Corollary 2 For any given connected graph G and every its induced subgraph G±, one has 
X m ax{G) > X max (Gi), so the synchronizability index of G satisfies r(G) < x mm (cl) ' ^ b°th G and 
G c are connected, then r(G) < , dmin . 

Since subgraphs have less nodes, this corollary is useful when a graph G contains some canonical 
graphs whose largest eigenvalues can be easily obtained as subgraphs (see the section below for 
further discussion). 

Corollary 3 For a given graph G, if the largest eigenvalue of G c is X max = d max (G c ) + a, 
then A 2 (G) = d m i n (G) + 1 — a. Consequently, the synchronizability index of G satisfies r(G) = 

dmin(G)+l—a ^ d mi „(G)+l-a 

^max(G) dmax (Ct) + 1 

By Lemma 2, generally a > 1, so the bound in Corollary 3 is better than the one in Corollary 2. 

3 Graphs having cycles as subgraphs 

In this section, an even cycle means a cycle with an even number of nodes and an odd cycle means 
a cycle with an odd number of nodes. By the discussion in the above section, one can get the 
following results. 

Lemma 5 For any even cycle, —2 is an eigenvalue of its adjacency matrix. For any odd cycle 

. / 3(ni — l)7r s 

sm( 2 x k — — ) 

with n\ nodes, — H (n "\w is an eigenvalue of its adjacency matrix. 

sin ( kr ) 



Proof Lemma 3 directly leads to the result. □ 
Theorem 1 For any given graph G, suppose G\ is its induced subgraph composing of all nodes 
of G with the maximum degree d max (G), and Q\ is the induced subgraph of G c composing of all 
nodes of G c with the maximum degree d max {G c ). Then, the following statements hold: 

(i) If both G\ and Q\ have even cycles as induced subgraphs, then X max (G) > d max (G) + 2 
and Xmax(G c ) > d max (G c ) + 2. Consequently, the synchronizability index of G satisfies r(G) < 

dmin(G) — l 
dmax(G)-\-2 

(ii) If both G\ and Q\ have odd cycles as induced subgraphs, and if the longest odd cycle of G\ 

Gin( 3(ni ~ 1)7r ) 

has n\ nodes and the longest odd cycle of Q\ has n-z nodes, then X max {G) > d max (G) + l ( ra 2 -iw 

sin(^ — — ) 

sm(-^| — J 

and \ m ax{G c ) > d max {G c ) + l ( " 2 1)7r . Consequently, the synchronizability index of G satisfies 

sinfi-* — — ) 



sin( ^ n — — ) 
dmin (G)H . („,_lW 

(iii) If Gi has an even cycle as an induced subgraph, and Q\ has odd cycles as induced sub- 
graphs with the longest odd cycle having n<i nodes, then X max (G) > d max {G) + 2 and X max (G c ) > 

sm(^-| — ) 

dmax{G c ) + 1 / "lw • Consequently, the synchronizability index of G satisfies 

V 2n 2 I 

. , 3(no —1)77 . 
sm( v g ' ) 

dmin(G)+ l (n 2 -Vj 

r(G)<- 



dmax (G) + 2 

(iv) If G\ has odd cycles as induced subgraphs with the longest odd cycle having n\ nodes 

s . n( 3(7^_l> ) 

and Q\ has an even cycle as an induced subgraph, then \ ma x{G) > d max (G) + 1 (n "\ )7r and 

sin(-^-75 — — ) 

v 2nj ' 

^max(G c ) > d max (G c ) + 2. Consequently, the synchronizability index of G satisfies 
sm < — hr, — > 



d max {G)+l — . (n 

Proof (i) By Lemmas 3 and 5, for any even cycle, its largest eigenvalue is 4 and -2 is an 
eigenvalue of its adjacency matrix. Let L\ be the sub-matrix of the Laplacian matrix of G related 
to all the nodes in G\, and A(G\) be the adjacency matrix of G\. Then, A{G\) contains a sub- 
matrix which is the adjacency matrix of the corresponding even cycle. Therefore, one has 

(d max (G) + 2)1 - Li = 21 + A(Gi) ? 0. 

In fact, the sub-matrix of 21 + A{G\) corresponding to the even cycle has a zero eigenvalue. This 
implies that the largest eigenvalue of G\ is larger than or equal to d max + 2. Then, by Lemma 1, 
Kiax(G) > d max + 2. Similarly, X max {G c ) > d max + 2. 

Further, suppose the number of nodes of G is N. Then, the smallest nonzero eigenvalue of G 
is A2 = N — \ ma x{G c ) and the minimum node degree of G is d m i n (G) = N — 1 — d max (G c ). If 
Xmax(G c ) > d max {G c ) + 2, then A 2 < N - d max {G c ) -2 = d min (G) - 1. Therefore, Theorem 1 (i) 
obviously holds. 



(ii) Follows easily from (i) and Lemma 5. 

(iii) and (iv) can be similarly proved. □ 
Remark 1 For a given graph G, clearly, the bounds given in Theorem 1 for the synchroniz- 

ability index r(G) are better than the general bound r(G) < ^ mi " . And the existence of cycles 
can be easily tested by drawing graphs. Since the networks with good synchronizabilities always 
have homogeneous degree distributions [10] , Theorem 1 is especially useful for estimating the 
synchronizabilities of homogeneous networks. 

Remark 2 The smallest even cycle is cycle C4, and its complementary graph C| has two 
separated edges. C4 and C| are very important in graph theory [15]. A graph has a C4 as an 
induced subgraph if and only if G c has C| as an induced subgraph. 

Remark 3 The cycle C3 is the smallest cycle. If a graph G contains only C3 cycles, its largest 
eigenvalue satisfies \ ma x{G) > d max + 1. 

Example 1 Consider cycle C5 in Fig. 1. Its complementary graph is also C5 (see Fig. 2). 
Testing the eigenvalues of C5 and its synchronizability, one finds that they attain the exact bounds 
in Theorem 1 (ii). 





Fig. 1 Graph C 5 Fig. 2 Graph C% 

Example 2 Consider graph Y\ in Fig. 3. Its complementary graph is C% (see Fig. 4). Testing 
the eigenvalues of Y\ and its synchronizability, one finds that they attain the exact bounds in 
Theorem 1 (i). 





Fig. 3 Graph T x 



Fig. 4 Graph = C 6 



4 Applications of other subgraphs 

According to the discussion in the above section, cycles are very important for estimating the 
network synchronizability. In this section, consider the estimation of the largest eigenvalue of a 
given graph by its subgraphs. First, the following lemma for adjacency matrices is needed. 

Lemma 6 [T7j Given a graph G, the smallest eigenvalue of its adjacency matrix A(G) satisfies 

^mm(^(G)) < d max (G) — X max 

By this lemma, one can get the following result. 

Theorem 2 For a given graph G, let H be a subgraph of G containing all nodes of G with 
the same node degree d. Let H\ be a subgraph of H. Then, the largest eigenvalue of G satisfies 

Xmax (G) > d + Xmax {H\ ) — d max (H\ ) . 

Proof It follows from Lemma 6 and the proof of Theorem 1. □ 

Remark 4 In the above section, cycles as subgraphs have the same node degrees. So the 
smallest eigenvalues of the adjacency matrices of cycles can be exactly computed. In Theorem 2, 
for a general subgraph Hi, its node degrees are not necessarily the same. In this case, Lemma 6 is 
very useful for estimating the largest eigenvalues. 

In what follows, consider some canonical subgraphs which can be used to estimate the largest 
eigenvalues of the original graphs. All the results can be combined with the results in the above 
section to obtain upper bounds for the network synchronizability. 

4.1 Graphs having chains as subgraphs 

A chain is a nearest-neighbor coupled graph without an edge between the first and the last nodes. 
For a given chain Pn with N nodes, it is well known that its largest and smallest nonzero eigenvalues 
are 2(1 + cos(-^)) and 2(1 — cos(-^)), respectively jHUHllII]- Obviously, the largest eigenvalue of 
a chain converges to 4 as its length tends to +oo. According to the above discussion, one can 
establish the following result. 

Theorem 3 For a given graph G, let G\ be a subgraph of G containing all nodes of G with the 
maximum degree. Suppose G\ contains a chain P^ with k nodes as its subgraph. Then, the largest 
eigenvalue of G satisfies X m ax{G) > d ma x + 2cos(^). 

Proof Theorems 1 and 2 directly lead to the conclusion. □ 

Note that chain P4 is a very important graph in graph theory [15]. The complementary graph 
of P4 is still a P4 graph. Therefore, a graph G has P4 as a subgraph if and only if G c has P4 as a 
subgraph. 



4.2 Graphs having bipartite graphs as subgraphs 

Let Gi(Vi, £\) and G2CV2, £"2) be two graphs on disjoint sets of n\ and nodes, respectively, where 
Vi and Si are the corresponding sets of nodes and edges. Their disjoint union G\ + G2 is the graph 
G\ + G2 = (Vi |J V2, £\ |J £%)■ Then, a bipartite graph G\ * G2 generated by G\ and G2 is a new 
graph composing of G\ + G2 and the new edges connecting each node of G\ to every node of G2, 
as detailed in [H [T5l [20] . It is well known that the largest eigenvalue of the bipartite graph G\ * G2 
is Knax{G\ * G2) = n\ + n2, i.e., the sum of the numbers of nodes of G\ and G2. In fact, the 
complementary graphs of bipartite graphs are always disconnected. 

By Lemma 6, for a bipartite graph G\ * G2, as discussed above, one can has \ m in(A(Gi * G2)) < 
dmax{Gi * G 2 ) - ni - n 2 . 

Theorem 4 For a given graph G, let i? be a subgraph of G containing all nodes of G with 
the same degree d. Suppose H contains a bipartite subgraph H\ * H2, and the numbers of nodes 
of Hi and H2 are ri\ and ri2, respectively. Then, the largest eigenvalue of G satisfies A max (G) > 
d + rai + n 2 - (i maa; (ill * if 2 )- 

Proof This is a direct result of Theorem 2 for bipartite subgraphs. □ 

Example 3 Consider graph T2 in Fig. 5. Obviously, T2 has a bipartite graph H as its subgraph 
which is composed of all nodes with degree 6 from T2- The largest eigenvalue of this bipartite 
graph is 8. So, by Theorem 4, A max (r2) > 9. On the other hand, the maximum degree of the 
complementary graph of T2 is 8. And the nodes with degree 8 in form a cycle G4. By Theorem 
1, the smallest nonzero eigenvalue of T2 satisfies A2(r2) < d m i n (T2) — 1 = 2. So, r(r2) < |. Simply 
computing the Laplacian eigenvalues of T2, one has A2 = 1.7251 and \ ma x = 9.2749. Consequently, 
r(r2) = q 2749 ~ 0.176. Therefore, the theorems in this paper successfully give the upper integer of 
the largest eigenvalue and the lower integer of the smallest nonzero eigenvalue of T2- 



1 5 




Fig. 5 Graph T2 



4.3 Graphs having product graphs as subgraphs 

Let Gi(Vi,fi) and ^2(1^2,^2) be two nonempty graphs. Their product graph G\ x G2 is a graph 
with node set V\ x V2, and (x\, X2)(yi, yi) is an edge in £{G\ x G 2 ) if and only if either x 2 = 2/2 
with x\y\ G £{G\) or xi = yi with X22/2 € £(G 2 ). One may view Gi x G2 as the graph obtained 
from G\ by replacing each of its nodes with a copy of G2 and each of its edges with the number of 
nodes of G2 edges joining the corresponding nodes of G2 in two copies, as detailed in [Tj [15]. 
It is well known that the largest eigenvalue of the product graph G\ x G2 is X max {G\ x G2) = 

X ma x{Gl) + ^max{G2)- 

Theorem 5 For a given graph G, let H be a subgraph of G containing all nodes of G with the 
same node degree d. Suppose H contains a product graph H\ x H2 as its subgraph. Then, the 
largest eigenvalue of G satisfies X m ax(G) >d + X max (Hi) + X m ax(H2) - d max (Hi x H 2 ). 

Example 4 Consider graph in Fig. 6. Obviously, all nodes of T3 have degree 4. And 
T3 has a product graph C4 x P3 (nodes 1 to 12) as its subgraph, where P3 denotes a chain with 
three nodes. The largest eigenvalue of this product subgraph is 7. So, by Theorem 5 or Lemma 
1) X max (r3) > 7. On the other hand, the complementary graph Tg of T3 has a bipartite graph 
as its subgraph which is composed of nodes 1 to 4 and nodes 9 to 12. By Theorem 4, similarly 
to Example 3, the largest eigenvalue of Tg satisfies A(r§) > d max (T^) + 3. Thus, by Corollary 3, 
the smallest nonzero eigenvalue of T3 satisfies A2(r3) < d m i n (T2) — 2 = 2. So, r(Ts) < |. Simply 
computing the Laplacian eigenvalues of T3, one has A2 = 1.2679 and A max = 7.4142. Consequently, 
r(T2) = 74142 ~ 0.171. Similarly to Example 3, the corresponding theorems proved in this paper 
successfully give the upper integer of the largest eigenvalue and the lower integer of the smallest 
nonzero eigenvalue of T3 . 
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Fig. 6 Graph T3 



4.4 The maximum disconnected subgraph 



Given a graph G of size N, suppose G\ is an induced subgraph of G of size n\, and G\ is dis- 
connected. G\ is called a maximum disconnected subgraph, if the node number of any other 
disconnected subgraph of G is less than or equal to n\ . 

Theorem 5 For a given connected graph G of size N, if the node number of its maximum 
disconnected subgraph is m, then the smallest nonzero eigenvalue of G satisfies A2 < N — n\. 
Consequently, r(G) < d Z~(G) + r 

Proof Suppose G\ is a maximum disconnected subgraph of size n\. Then, by Lemma 4, 
\nax{G\ ) > m. Further, X m ax{G c ) > n 1} so A 2 (G) < N - m. □ 

Example 5 Consider graph T4 in Fig. 7. By deleting node 3 or 6, one can verify that the node 
number of its maximum disconnected subgraph is 7. So, A 2 (r4) < 1. Combining Lemma 2 and 
Theorem 6, one has r(T^) < |. It is well known that graphs as in Fig. 7 have large node and edge 
betweenness centralities therefore have bad synchronizabilities. Theorem 6, based on the theory 
of subgraphs and complementary graphs, gives an explanation why such graphs indeed have bad 
synchronizabilities . 




Fig. 7 Graph T 4 

Remark 5 From Examples 3 to 5, it can be seen that actually one does not need to draw 
out the complementary graphs; the properties of complementary graphs can be analyzed using the 
original graphs. 

5 Conclusion 

In this paper, some relatively tight upper bounds have been obtained for the network synchroniz- 
ability based on the theory of subgraphs and complementary graphs. Especially, some sharp bounds 
are given for the eigenratioes and the largest eigenvalues when a graph has some canonical sub- 
graphs, such as cycles, chains, bipartite graphs and product graphs. Considering that the networks 
with good synchronizabilities typically have homogeneous degree distributions, the results obtained 
in this paper are particularly useful for estimating the synchronizabilities of homogeneous networks. 



And for a given network, its corresponding graph as well as its complementary graph and subgraphs 
can be easily obtained. This paper shows that better understanding and helpful manipulation of 
subgraphs and complementary graphs are very useful for enhancing the network synchronizability. 
Therefore, the graph-theoretical method provided in this paper is deemed important in the study 
of network synchronization problems. 
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